Prom microscopic to macroscopic dynamics in mean-field theory: effect of neutron 

skin on fusion barrier and dissipation. 
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In this work, we introduce a new method to reduce the microscopic mean-field theory to a classical 
macroscopic dynamics during the initial stage of fusion reactions. We show that TDHF (Time- 
dependent Hartree-Fock) could be a useful tool to gain information on fusion barriers as well as 
on one-body dissipation effects. We apply the mean-field theory to the case of head-on reaction 
between ie O and 16 > 22 > 24 - 28 o in order to quantify the effect of neutron skin on fusion. We show 
that the determination of fusion barrier requires, in addition to a precise knowledge of the rel- 
ative distance between the center of mass of the two fusing nuclei, the introduction of an addi- 
tional collective coordinate that explicitly breaks the neutron-proton symmetry. In this context, 
we estimate the position, height and diffuseness of the barrier as well as the one-body friction 
and show that a global enhancement of the fusion cross-section is expected in neutron rich nuclei. 

(February 6, 2008) 
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I. INTRODUCTION 

The understanding of the formation of compact sys- 
tem during fusion of two nuclei implies the development 
of nuclear models which can take into account quantum 
effects, the dynamical evolution and dissipation in a com- 
mon framework. From that point of view, the description 
of the conversion of relative motion into internal excita- 
tion of the composite system and also the deexcitation of 
the compound nucleus is far from being complete. Due to 
the diversity of dissipative mechanisms, the microscopic 
descriptions of fusion reactions is often substituted by a 
macroscopic dynamics in which only few relevant degrees 
of freedom are considered. Such macroscopic description 
are only valid under a series of assumption like separa- 
tion in time-scales between relevant and other degrees of 
freedom. The validity of these hypothesis as well as the 
link with microscopic theories is under debate. 

In the present work, starting from quantum mean-field 
theory, we address the problem of the reduction of the 
microscopic dynamics to a macroscopic theory Jl]-|j|. 

A. Mean-field theories 

We consider the Time-Dependent Hartree-Fock (TDHF) 
theory: one-body degrees of freedom are supposed to be 
dominant and the description of the system reduces to 
the time evolution of a Slater-determinant \*f?) in a self- 
consistent mean-field. Noting p the associated one-body 
density matrix, the evolution is given by the Liouville 
equation M 



(1) 



tions, and gives reasonable results in the case of fusion 
reactions j7|-[)]]. ^ n addition, it has the advantage to take 
into account shell effects and to include the so-called one- 
body dissipation. 

However, fusion reactions and more generally deep in- 
elastic reactions are often described assuming that some 
macroscopic variables are dominant Although a 

large number of work is devoted to this subject, the 
link between macroscopic dynamics and microscopic ap- 
proaches such as TDHF is still unclear [jll]] . 

In this paper, we address this problem considering 
head-on reactions. In this case, the relative distance R 
between nuclei in the entrance channel is known to de- 
scribe the main features of the dynamics and its evolution 
follows the classical Hamilton- like equation Q || : 



dp 

di 



F(R) -j(R)R 



(2) 



where fJ>(R) is the reduced mass of the system and the 
force is assumed to derive from a potential ( F(R) = 
— The dissipation kernel j(R) is supposed to re- 

produce the effect of the other non described degrees of 
freedom on the considered variables. 

The dissipative effects are generally taken by including 
the one-body dissipation in a liquid drop picture |l3| or 
using an adiabatic approximation in a two-center shell 
model (see for instance [[l4|). However, the understand- 
ing of dissipative effects in nuclear collisions is far from 



where h [p] is the self-consistent field. Such a picture has 
been widely used to investigate nuclear dissipative reac- 



*For simplicity, we have written the expression in its sim- 
plest form, it should however be noted that other degrees of 
freedom like angular momentum ji or mass assymetry ]l^ ] 
may be considered in the formalism. 
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being complete and the use of equation (g) rises some 
questions. A large number of effects are known to par- 
ticipate to the dissipation of the relative to the inter- 
nal excitation energy during the reaction. Among them, 
exchange of particles and deformation are generally as- 
sumed to be dominant and are accounted as one-body 
dissipation effects. We would like however to mention 
that the dissipative process described by equation (Q) is 
obtained using quantum perturbation theory Q. It is 
not clear whether it can be applied to the case of violent 
collisions like fusion. 

Mean-field theory like TDHF does a priori includes 
one-body dissipation in a non-perturbative framework. 
In addition, it assume neither separation between degrees 
of freedom nor adiabaticity. However, the link with the 
so-called one-body dissipation theory is still under dis- 
cussion fl6|] . 

It should also be noted that friction models assumes 
a fast equilibration of other degrees of freedom as com- 
pared to the relative motion leading to the formation of a 
compound nucleus at finite temperature. Due to the ab- 
sence of two-body or higher order correlations, TDHF is 
unable to account for this effect. From the experimental 
point of view, the equilibration time is estimated to be 
of the order of the reaction time |l7]] , which means that 
the classical description assuming thermal equilibration 
assumption is questionable as far as fusion is concerned. 

In this paper, we discuss the possibility to use macro- 
scopic equations for TDHF although it does not assume 
perturbation theory or equilibration. In the following, 
we show that such a reduction seems possible and gives 
reasonable results on fusion barriers properties. 

Different works already exist on the possibility of such 
a reduction of microscopic theory fll8 25-2^|. However, 
as it was already noted in p8[ , high numerical accuracy 
as well as the use of complete realistic effective forces 
is required. In this article, we use the three-dimensional 
TDHF code developed by Bonche and co-workers [|lQf| and 
the full SLy4d Skyrme force including spin-orbit pcf . 

In the following, we apply the mean-field theory to 
head-on reactions between oxygen isotopes. The static 
properties of oxygen isotopes are first presented. We 
then describe the method used to reduce mean-field evo- 
lution to the relative distance only and extract quantita- 
tive information on the fusion process. The macroscopic 
potential landscape and the dissipative aspects are then 
discussed. 



II. RESULTS 

A. Structure of oxygen isotopes 

For each of the considered nuclei, the self-consistent 
Hartree-Fock equation reads: 



h[p] yi) = e % \<pT) 



(3) 



where h[p] is the mean- field hamiltonian, | ipj) denotes the 
single particle wave functions and r denotes the isospin 
quantum number (r = p [proton] or r = n [neutron]). 
Equation [|is solved in r-space on a mesh of size (16/m) 3 
with a step of 0.8 fm. The complete expression of h[p] 
for Skyrme forces can be found for instance in plfl . The 
resulting single-particle levels energies obtained using the 
SLy4d Skyrme force are presented in Fig. [l] for the four 
isotopes 16 0, 22 0, 24 and 28 0. 
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FIG. 1. Proton (p) and neutrons (n) energy levels obtained 
in Hartree-Fock calculation described in the text using the 
Sly4d effective Skyrme force for 16 O, 22 O, 24 O, 28 O isotopes. 

The associated neutron and proton root mean square 
radius (RMS) are reported in figure |^. These root mean 
square radii leads to RMS matter radii in agreement with 
experimental data for the 16 0, 22 and 24 nuclei p8[ . 
A neutron skin is evidenced when the number of neutrons 
increases. 
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FIG. 2. Proton (circle) and neutron (diamond) root mean 
square radius predicted by the Hartree-Fock calculation us- 
ing the Sly4d Skyrme interaction for the 16 0, 22 0, 24 0, 28 
isotopes. 

Once the static properties are obtained, we consider 
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the evolution of the two colliding nuclei in a larger vol- 
ume (32/m) 3 with the same step in r-space as before. 
The initial distance between the center of the two nuclei 
is equal to 16 fm. The time step is chosen as 0.45 fm/c. 
Before showing the results of the dynamical calculation, 
we discuss the method used to reduce the mean-field dy- 
namics to the time evolution of a few macroscopic coor- 
dinates. 



B. Reduction to macroscopic dynamics 

In order to define the relative distance between the 
two nuclei, we consider the mean-field equation in r-space 
and the associated one-body density p T (r,r') with p T — 
E» 1^1 W) I- We assume, as in rcf. Q, that at each 

time step, the system could be split into two subspaces 
1 and 2 with the associated densities p r x (x = 1,2). We 
define the separation plane after the contact at the neck 
pos ition as it is done generally in macroscopic models 

Once the separation is obtained, we compute for each 
subspace, the number of neutrons N x — Tr (p™) or pro- 
tons Z x = Tr (p™), the mean position 



Rl=<r> x =Tr {fp r x )/Tr{pl) 



and the associated momentum 



Pl=<p> x =Tr{ppl) 



(4) 



(5) 



where Tr denotes the Trace. Let us first assume that we 
can neglect the isospin degree of freedom (the validity of 
this assumption will be discussed in the following). We 
define the center of mass as well as the associated total 
momentum: 



R x = (N X R X + Z X RP X ) /A x 
P x = PL 1 + PS. 



(6) 
(7) 



where we have introduced A x = N x + Z x . 

Using Hamilton equation for each subsystem, we estimate 

the mass from equation 



PJ 



dR x 
dt 



(8) 



At a given relative distance R = R\ — R 2 , we compute 
the reduced mass associated with the composite system 
p(R) = mi 7712/ (mi + 7712). Note that, since we have 
a Skyrme force, with in particular t\, t 2 and spin-orbit 
terms, p(R) is different from the normal value: 



MR) 



A 1 {R)A 2 (R) 
(Ai(R) + A 2 {R)) 



(9) 



where m is the nucleon mass. In our calculation, wc 
obtain p(R) / po{R) — cte ~ 0.67 for all systems. This 
remains valid up to the formation of a compact shape. 





Once these quantities have been defined and calcu- 
lated, the different components entering in Eq. (||): R(t), 
dR(t) /dt, P{t) and its derivative can be computed. 

Eq. ^contains two unknown quantities, the force F(R) 
and the friction coefficient 7(i?). In order to estimate 
these quantities, a method proposed in ref. |24l , [l8|] is used. 
It consists in performing mean-field evolution with two 
slightly different initial energies E 1 , E 11 . At a given 
position R, Eq. || can be interpreted as a system of two 
equations with two unknown quantities which reads: 



(10) 



where all derivatives are taken at position R. We checked 
the numerical stability of the proposed procedure, by 
solving the classical equation associated with the same 
initial condition as TDHF evolution and with the fric- 
tion and potential extracted from it. As expected, we 
recover the TDHF evolution. 



C. Nuclear potential 

We performed TDHF calculations with four different 
initial center of mass energies E/A = 0.60, 0.65, 0.70 and 
0.75 MeV for each reaction. Figure || gives the poten- 
tial for the reactions 16 0+ 16 ' 22 ' 24 ' 28 obtained by inver- 
sion of Eq. [ITj. To extract the nuclear part, we assume 
V nuc i ear (R) = V(R) — Vc{R), where the coulomb part is 
simply taken as V C (R) = e 2 Z 1 (R)Z 2 (R)/R (dashed line). 
For each reaction, we have solved Eq. 10 using different 



couple of TDHF evolution with different initial energies. 
The solution depends slightly on the couple of trajecto- 
ries chosen, this uncertainty is represented by errorbars 
in Fig. H The uncertainties are rather small indicating 
that the microscopic dynamics is already well described 
byEq. I 

Examining the global shape of V(R), and in particular 
the maximum, we see that although the barrier height is 
slightly affected by the neutron skin, the increase of the 
neutron number not only increases the barrier position 
Rb but also change its diffuseness. This effect should 
normally indicate a very different behavior in the sub- 
barrier fusion of very neutron rich nuclei as compared 
to stable isotopes. Using a different technique, the same 
conclusion is drawn in ref. [^9| while in ref. (3^Jl^], it 
is concluded that the excess of neutron does not affect 
fusion. This point is discussed in the following 
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In order to compare different systems, "f{R) and 
Ei oss (R(t)) are plotted as a function of the relative dis- 

r\ > x l 2 + <r 2 p 



tance normalized to R\2 = I (< r% > 1 / 2 + < r% > 
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FIG. 3. Potential land- 

scape obtained in head-on ie O+ 16,22 ' 24 ' 28 reactions. The 
total (thick solid lines), coulomb (dashed line) and nuclear 
potential (thin solid lines) are presented. For the calcula- 
tion, different couples of trajectories taken within the initial 
energies E/A = 0.75, 0.70, 0.65, 0.60 MeV are presented. "Er- 
rorbars" correspond to the mean-variance on results obtained 
with different initial conditions. The isodensity plots in the 
reaction plane are also superimposed for different times. 



D. Dissipation 

While numerous works have been devoted to poten- 
tials determination, only few quantitative analyses exist 
on the dissipation effects starting from a non-adiabatic 
microscopic theory. The method described above gives 
also information on the dissipative processes. The fric- 
tion coefficient associated to the mean-field evolution of 
the system can be computed. It is presented in Fig. |J. 

In order to have a global information on the amount of 
energy transferred from the relative motion into internal 
degrees of freedom of the compound system, the inte- 
grated value of the Rayleigh function can be evaluated. 
The Rayleigh function gives the rate of energy dissipated 
by unit of time: 



E ray (t) = dE/dt = -~f(R)R 2 



(11) 



and its integrated value Ei oss (R(t)) reads (Fig. [|): 

E loss (t) = 2 f E ray (s)ds (12) 

Jtn 



where P and T refers respectively to the projectile and 
target. It is worth noticing that dissipation is larger 
at larger distances than R/R12 ~ 0.7 in the case of 
neutron rich nuclei. However, at intermediate distance 
0.5 < R/R12 < 0.7 the dissipation is larger in the sym- 
metric case. We indeed expect a different behavior be- 
tween symmetric and assy metric systems. If the one- 
body dissipation picture is correct, part of the viscosity 
is due to the change of the global shape of the system (the 
so-called wall dissipation 31|). In assymetric reactions, 
this effect should be completed by the friction due to the 
exchange of nucleons between the two subsystems (the 
window dissipation p^| ). The link between the friction 
coefficient in TDHF and the wall plus window picture 
is however not clear. Indeed, mean-field theory also in- 
cludes other degrees of freedom which are not taken into 
account in the wall and window formula. In particular, a 
fraction of the initial energy is converted into collective 
excitation (giant resonances). In our case, this effect is 
also integrated in "f(R). 

We finally would like to mention that, although the 
dissipation is different between symmetric and assymet- 
ric systems, the value of Ei oss at minimal distance of 
approach is independent of the assymetry. 
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FIG. 4. Top: Friction coefficient as a function of relative 
distance of head-on 16 0+ ie O (circles), 16 0+ 22 (crosses) and 
16q_|_28q ^hick lines) reactions. Bottom: Associated inte- 
grated Rayleigh function Ei oss as a function of the relative 
distance normalized to the sum of the nuclear radius of the 
target and projectile R12 = \/f (< »4 > 1/2 +<r% > 1/2 ). 
In these figures errorbars are not shown. For simplicity, we 
do not show the 16 0+ 24 which is similar to the 16 0+ 22 



E. Precise determination of the barrier properties 

The possibility to obtain fusion cross-section from 
mean-field requires a precise estimate of the nuclear po- 
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tential in the vicinity of the fusion barrier. Details of the 
fusion barrier extracted using Eq. || are shown in figure 
H (gray area). 
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FIG. 5. Zoom of the fusion barrier region obtained for 
16 0+ 16,22,24 ' 28 reactions respectively from top to bottom us- 
ing the macroscopic interpretation of TDHF. Filled areas cor- 
responds to calculated potential with uncertainties, the gray 
area corresponds to the result obtained using Eq. H while the 
black area is the result obtained including the isospin relative 
distance difference (see text). In the case of 16 0+ 16 there 
is a simple curve since the additional degrees of freedom does 
not change the extracted potential. 

For 16 0+ 22 and 16 0+ 24 reactions, rather large un- 
certainties are observed for specific relative distances. 
These uncertainties represents the dependence of the po- 
tential V(R) with the initial energies, i.e. the polarisa- 
tion of the extracted solution using Eq. ^ with mean-field 
trajectories. In our microscopic calculations, this polar- 
isation could be directly attributed to the appearance 
of another important degree of freedom which is not in- 
cluded in Eq. [| Indeed, regions where uncertainties are 
important are associated with relative distances between 
protons in the two nuclei differing from the distances be- 
tween neutrons. In order to account for this additional 
effect, we have introduced the difference between the pro- 
ton and neutron relative distances: 



AR np = - R») - {R{ - Rl) (13) 

and assumed that the Eq. |^ is modified in order to in- 
clude this additional degree of freedom: 

dP 

— = F(R)- 7 (R)R- lnp (R)AR np (14) 

As in the previous case, using three different initial condi- 
tions, Eq. [li] can be interpreted as three linear equations 
with thre e u nknown quantities F(R), 7 (R) and 7„ p (i?) 
(see Eq. [Ad]). The inversion of such ensemble of equa- 
tions is rather difficult since the additional degree of free- 
dom is not always important during the evolution. This 
often leads to an overcomplete set of equations. A dis- 
cussion on inversion method is given in appendix A. 

Assuming that the dynamical mean-field theory can be 
reduced to Eq. |lj, we obtain the new potentials repre- 
sented by a black area in Fig. g. In that case, the initial 
energy dependence are strongly reduced as compared to 
the previous calculation indicating the importance of this 
new degree of freedom. Note that in the symmetric re- 
action 16 0+ 16 0, the additional degree of freedom can 
be neglected, and the obtained potential strictly corre- 
sponds to the old one. It should be noticed that the 
second maximum present in 16 0+ 22 and 16 0+ 24 re- 
actions using Eq. |^ has disappeared. In order to quantify 
more precisely the effect of the neutron excess on fusion, 
we have fitted the new value of V(R) in the vicinity of 
the barrier assuming V{R) =V B - AV B (R - Rb) 2 - We 
have reported in table |, the barrier height Vb , the peak 
position R b and th e deduced curvature Hu>b defined as 
huj B = h^AV B /na- 

In the case of 16 O + 16 O, where experimental data ex- 
ist, a good agreement with the extracted barrier height 
and position f p3p^ ] is found. For other reactions, ex- 
perimental results are not available. We would like to 
mention that we have also applied this method to the 
formation of heavy and very heavy elements. In all cases, 
we found a good matching with the experimental fusion 
properties giving additional confidence in the physical 
meaning of the extracted barriers. 

Looking at the three assymetric reactions, the barrier 
height is almost constant while, as expected, the peak 
position increases slowly as well as the barrier diffusion. 
However these trends could not be extended to the sym- 
metric case. This fact associated to the fact that it is not 
necessary to introduce additional degrees of freedom in 
the symmetric reaction case, shows the peculiarity of the 
neutron excess in fusion reactions. 
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24 O 
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Vb (MeV) 


9.8 


8.9 


8.8 


8.9 


Rb (fm) 


8.6 


9.6 


9.8 


10.0 


hulB 


2.39 


1.21 


1.36 


1.65 



TABLE I. Fusion barrier Vb and position extracted from 
head on collisions i6o+ 16 ' 22 ' 24 ' 28 0. 
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In order to have some qualitative information on fusion 
cross sections, we have injected the values reported in 
table I in the Wong formula : 



a F (E) 



huj B R 2 B 
2E 



In (1 + exp (2ir (E - V B ) /nu B )) (15) 



where E is the laboratory energy. The cross sections for 
16 0+ 16q ( thin Une ) j i6 0+ 22 (circles) and 16 0+ 28 
(thick line) reactions are displayed in Fig. ^. For fu- 
sion above the barrier, due to the 1 MeV lower value 
of the barrier, a larger fusion cross section for the case 
of neutron rich isotopes as compared to symmetric 16 O 
reaction is expected. However, there is no differences 
between 22 0, 24 or 28 above E > lOMeV. 

The excess of neutron is also predicted to modify the 
sub-barrier cross-section [^0|,|6[|7): an increase of the 
sub-barrier cross-section associated with a decrease of 
the slope of the cross-section is expected. In Fig. ||, 
a more complex dependence in the slope of the cross- 
section is observed. It should however be noted that the 
inclusion of friction, which is neglected in the Wong for- 
mula, may change the slope in the sub-barrier region. 
Since the friction changes considerably between stable 
and neutron rich nuclei (see Fig. |j) , a comparison of the 
different slopes is difficult. In the case of the two reac- 
tions le O + 22 O and le O + 28 O, friction coefficients have 
the same behavior as a function of the relative distance 
and the cross-sections may be compared. In this case, we 
indeed observe a smaller slope in the 28 case. 
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FIG. 6. Example of calculated fusion cross-section using 
the Wong formula for the reaction 16 0+ 16 (thin line), 
ie 0+ 22 ( circles ) and i6o+ 28 (thick line). 

The evolution of 7„ p (i?) as a function of the relative 
distance is shown in Fig. [?] as well as the position of the 
barrier. The effect of the new coordinate is intermittent. 
The first maximum of 7„ p is close to the barrier. This 
is directly understood by the fact that, at this point, 
the difference between proton and neutron distances is 
maximum. Indeed, in this case, the coulomb repulsion 
is large, which goes in favour of a larger proton distance 
while, due to the existence of a neutron skin, the attrac- 
tive nuclear part acts first on neutrons and thus tends 



to reduce the neutron distance. After the contact, the 
effect of 7 np appears periodically due to the appearance 
of neutron-proton oscillations which are already present 
in the entrance channel. 

The necessity to introduce explicitely the difference be- 
tween neutron and proton center of mass in neutron-rich 
nuclei has already been noted in several work [[39] [IT[ . 
In all cases, it was concluded that the neutron excess 
is expected to increase the fusion cross-section as also 
observed in the present work. Different effects are in- 
deed expected to act in this direction: the large exten- 
sion of the neutron skin which implies the influence of 
the nuclear field at larger relative distance, the coupling 
between the soft dipole mode with the relative motion. 
The time-dependent mean-field model includes all these 
effects. In particular, the observable AR np is obviously 
connected to the polarisation of each nucleus before fu- 
sion. It should also be noted that AR np is affected by 
the transfer of particles initiated just after the contact in 
order to equilibrate the N/Z ratio between the two part- 
ners. This transfer which is known to induce damping in 
the relative motion, gives rise to the excitation of a pure 
collective GDR in the composite nucleus [Q . 

Although the additional friction term in Eq. (14) 
seems to strongly reduce uncertainties on the extracted 
potentials, for longer time, a more complex dynamics is 
expected. In particular, it would be convenient to treat 
also collective processes like the excitation of giant reso- 
nances. For instance, the effect of polarisation as well as 
the onset of the GDR after fusion requires the introduc- 
tion of potential effect coming from the symmetry energy. 




FIG. 7. Top: Friction coefficient 7„ p as a function of the 
relative distance in head-on 16 0+ 22 (top), 16 0+ 24 (mid- 
dle) and 16 0+ 28 (bottom) reactions. The barrier position 
is indicated by an arrow. 

To quantify the amount of energy transferred from the 
relative motion to the new coordinate, we have also com- 
puted a quantity similar to Ei oss noted E np which is de- 



G 



fined as: 

E np (R) 



lnp (R(s))R{s)AR np ( s )ds (16) 



The result of the computation is shown in Fig. ^. The 
amount of transferred energy to the new observable is 
directly proportional to the neutron excess. 



2.0 




FIG. 8. Energy transferred between the relative mo- 
tion and the neutron to proton oscillation as a function of 
the relative distance for the reaction 16 0+ 22 (dashed line), 
16 0+ 24 (thin line) and 1B 0+ 28 Q (thick line). 



III. CONCLUSION 

In this paper, we have shown that the properties of 
a macroscopic dynamics can be rather precisely inferred 
from a microscopic quantum mean-field dynamics by re- 
ducing the one-body information to few collective coor- 
dinates. In order to have quantitative information on 
the barrier properties, the "naive" reduction taking into 
account only the relative distance is not sufficient in neu- 
tron rich nuclei and it is necessary to introduce the dif- 
ference between neutrons and protons during the onset 
of fusion. 

We have applied our method to the case of head-on 
oxygen reactions involving oxygen isotopes. We have 
shown that, while the fusion barrier is only slightly af- 
fected by the neutron skin, the barrier shapes change. 
We thus expect that the neutron skin will modify prefer- 
entially the sub-barrier fusion reactions. 

Beside the potential landscapes, the mean-field dynam- 
ics gives also information on how the relative motion dis- 
sipates into internal excitation of the composite system. 
Such a dissipation which gives directly the associated fric- 
tion coefficient has rarely been obtained in non-adiabatic 
microscopic theory. We have shown that during the ini- 
tial stage of the reaction, part of the dissipation corre- 
sponds to an isospin symmetric dissipation of the relative 
motion (7). In this case, the total amount of energy de- 
posited in the system seems independent of the neutron 
skin. On the other side, another part of the dissipation 
associated with the neutron-proton difference (7« p ) de- 
pends explicitely on the neutron excess. 



In this work, we have shown how mean-field theory 
which includes nuclear structure effects and does not in- 
voke specific approximations on macroscopic dynamics, 
may be very useful to extract quantitative information 
on the potential landscapes as well as on dissipation. 

The techniques developed in this paper may be of par- 
ticular interest for super-heavy systems in which macro- 
scopic theories are presently widely used p8| , ^3| -^5|. 
There, friction coefficients are often calculated assuming 
thermal equilibrium of other degrees of freedom. Such 
an assumption is not included in TDHF and it would be 
desirable to compare the two approach. 
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APPENDIX A: INVERSION OF EQUATION 
WITH AR NP 

Using three TDHF trajectories, noted ( I, II and III), 



Eq. 14 can be interpreted as a matricial equation which 



reads for a given position R: 



1 - 



Vi - 



dR 1 


dAR' np 


dt 


dt 


dR 11 




dt 


dt 


dR 111 




dt 


dt 



7(i?) 

\ lnp(R) I 



I dt \ 



dP 1 



dt 



dP' 



(Al) 



dt 



which we write M.X = B in the following. The inversion 



of Al is made difficult due to the fact that the additional 



degree of freedom acts intermittently on the other com- 
ponents. When the coupling between AR np of freedom 
is small, the set of equations becomes overcomplete. 

In order to solve this difficulty, we first consider 
two TDHF trajectories and solve the two-dimensional 
Eq. [l0| This gives an approximative solution of Al 
(Fo(R), 70 (i?)). A measure of the importance of AR np 
is then obtained by computing the residue of the three 
dimensional equation: 



Res{R) = \\MX -B\ 



(A2) 



where X =* (F (R), 7o(-R), 0). Example of calculated 
residue for the 16 0+ 22 and 16 0+ 24 reactions are pre- 
sented in Fig. ^[ We see that a large residu at a given 
time (or equivalently position) initiate the errors for fu- 
ture time (or lower position). This comes from the fact 
that we compute V(R) by integration of F over R from 
large to small R values. 
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o 

X 



FIG. 9. Example of calculated residue for the reactions 
i6q + 22q ( tQp ^ i6q + 24q ( bottom ) as a f unct ion of the rela- 
tive distance. The potential obtained by inversion of equation 
is also superimposed with errorbars. 

In this figure, the appearance of larger and larger 
residues during the approaching phase towards the fused 
systems is clearly observed. This effect may be seen as 
the starting of neutron-proton oscillations. We also re- 
mark that in between the ma xima, the residue is very 
close to zero and equation Al reduces to Eq. [k| In our 
numerical imp leme ntation, when the residue is large, we 
solve equation Al while for very small values, we switch 
to Eq. p. 
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